Current oscillations in a metallic ring threaded by a time-dependent magnetic flux. 
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We study a mesoscopic metallic ring threaded by a magnetic flux which varies linearly in time 
$ M (t) = $t with a formalism based in Baym-Kadanoff-Keldysh non-equilibrium Green functions. 
We propose a method to calculate the Green functions in real space and we consider an experimental 
setup to investigate the dynamics of the ring by recourse to a transport experiment. This consists in 
a single lead connecting the ring to a particle reservoir. We show that different dynamical regimes 
are attained depending on the ratio h$/$oW, being $o = hc/e and W, the bandwidth of the ring. 
For moderate lengths of the ring, a stationary regime is achieved for ft$/$o > W. In the opposite 
case with h&/&o < W, the effect of Bloch oscillations driven by the induced electric field manifests 
itself in the transport properties of the system. In particular, we show that in this time-dependent 
regime a tunneling current oscillating in time with a period r = 27r$o/$ can be measured in the 
lead. We also analyze the resistive effect introduced by inelastic scattering due to the coupling to 
the external reservoir. 
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I. INTRODUCTION 

-a 

In the last years, the progress of technology has enabled the experimental study of very interesting phenomena 
taking place in mesoscopic systems. The fascinating feature in these systems is that several physical processes that 
have a pure quantum-mechanical origin can be observed in a suitable transport experiment. Paradigmatic examples 
are the Aharanov-Bohm effect in systems with annular geometry and Bloch oscillations in periodic structures. 
. The Aharanov-Bohm effect in circular devices motivated a vast theoretical and experimental activity for some time 
now. A metallic mesoscopic ring threaded by a static magnetic field displays persistent currents, which oscillate as a 
function of the magnetic flux <£>a/, with a period <I>o = 2iTh/e Like other mesoscopic phenomena, these persistent 
currents are a consequence of the coherence of the wave function throughout the system. The Aharanov-Bohm effect 
also manifests itself when such a ring is connected to two metallic leads at different external voltages. This geometry 
£S| ' is usually employed in experiments to measure the conductance of the metallic ring in the presence of the external 
flux. In this case, quantum interference cause oscillations in the current through the collector lead with a period <f>o 
as $m varies (|]. 

More than 70 years ago, Bloch predicted the occurrence of a time-periodic motion of the electrons in a band under 
the influence of an external constant electric field F j|. Recently, Bloch oscillations with a period h/eFL has been 
observed in superlattices being L the periodicity of the lattice and e the electron charge. This motivated new 
£^ ' theoretical development in the field ^-0]. 

A metallic ring penetrated by a time-dependent magnetic field is a very interesting system, worth of theoretical 
and experimental study. Since time is the unavoidable variable in any experimental measurement of currents while 
, the magnetic field is varied, that device is more realistic than its static-field counterpart. The case of a magnetic flux 
with a linear time dependence is particularly appealing because this situation is closely related with that of electrons 
in a periodic lattice under the influence of a dc voltage. In fact, a magnetic flux that varies linearly in time induces 
a constant electric field F = (1/ cL)d$> M(t) / dt, being L the longitude of the ring, while the geometry of the ring 
produces a periodic structure for the electrons inside it. While a good amount of literature can be found on the 
Aharanov-Bohm effect, considerably fewer studies have been reported on this problem. It has been addressed in an 
early work by Biittiker, Imry and Landauer ]8| who predicted that the electrons in an ideal one-dimensional ring 
should execute Bloch oscillations driven by the induced electric field. Some time later, Landauer and Biittiker showed 
that inelastic scattering in that system generates a dc current and dissipation qJ. A single lead connected between 
the ring and an external reservoir was later proposed by Biittiker as a concrete way to introduce inelastic scattering 
events |lC| ]. Since then, the attention focused mainly in Aharanov-Bohm loops where the flux threading the ring is 
static and devices with inserted quantum dots seem to occupy the center of the stage pdj|. 
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In this work the problem of the normal-metal loop in the presence of an external time dependent magnetic flux 
with the form = 3>t is studied in the framework of non-equilibrium Green functions. We consider an ideal ring 

described by a one-dimensional tight-binding model with a bandwidth W and we propose an approach to calculate 
the Green functions in real space. This procedure enables to gain insight on the different energy scales that intervene 
in the dynamics of the problem. In addition, it has the advantage of being also adequate to tackle more general 
situations where the configuration is not translational invariant, like the case of rings with disorder and with inserted 
barriers or quantum dots. 

The experimental study of the physical properties of the ring in the presence of the variable magnetic field requires 
some kind of coupling between this system and the external world. We consider a concrete experimental setup to 
study time-dependent charge oscillations in the ring. This consists in a single lead connecting it to a particle reservoir. 
The goal is to study the quantum tunneling between the ring and the reservoir, originated by the time dependent 
oscillations driven by the induced electric field. The metallic wire is described as an ideal system of non-interacting 
electrons and the reservoir is assumed to fix the chemical potential /i. We use an approach based in Baym-Kadanoff- 
Keldysh formalism to solve the coupled system. This method has been successful to derive generalized transport 
equations [jl2|,^3) in systems with an explicit time dependence. In the present context the main idea is to represent 
the effect of the lead and the reservoir by the self-energy effects they produce. 

In the case of pure one-dimensional rings pierced by a static magnetic field in the absence of disorder, the persistent 
currents vanish as 1/L, being L the circumference of the ring |E3|. In the time dependent case there is a non trivial 
connection between the size of the system and its dynamical behavior. We study the behavior of the retarded Green 
function as a function of $ and L. We argue that for a magnetic field varying at a rate such that h$/$o < W, the 
time dependent features of the relevant physical magnitudes vanish as 1/L. Instead, for more rapid changes of the 
magnetic field such that /i$/$o > W, their time dependence vanish exponentially fast as a function of L. Our aim 
is to show that for moderate lengths L and field variations satisfying the condition h<f>/<f>o < W, quantum tunneling 
takes place between the ring and the reservoir and a time dependent current could be measured in the lead. We also 
analyze the effect that inelastic scattering originated by the coupling to the reservoir cause on the transport properties 
of this system. 

The paper is organized as follows. Section II is devoted to expose the theoretical approach. In Section III we obtain 
the set of equations to calculate the retarded Green function in the ring unconnected from the reservoir as well as 
connected to it. Results are shown in Section IV. We finally summarize our study and present the conclusions in 
Section V. 



II. BASIC THEORY 



We consider an experimental setup to study the dynamics of the ring in the presence of the time-dependent magnetic 
flux of the form <S>m — $i. The device is sketched in Fig. la. Our aim is to study the quantum tunneling between the 
ring and an external particle reservoir by measuring the current flowing through the lead. We use an approach based in 
Baym-Kadanoff-Keldysh non-equilibrium Green functions to solve the coupled system. We follow a procedure similar 
to that used in Refs. [|l2|jl3| , where the effect of the lead and the reservoir is represented by a suitable self-energy. 
The system is simplified as indicated in Fig. lb. The lead and the reservoir are described as an effective system a 
with a fixed chemical potential fi. A contact, described by a single bond, establishes the connection between the ring 
and the lead. 

The total hamiltonian of this system contains three pieces: 

H = H r i n g + H a + H cont act- (1) 

The first term describes the metal ring threaded by the magnetic flux. We consider non-interacting spinless electrons 
in a tight-binding chain of longitude L with N sites, lattice parameter a = L/N and periodic boundary conditions 
(N + 1 = 1). The model hamiltonian is 

JV 

H ring = -T^Ce-^Vx + e^cJ^d), (2) 
/ 

with I — 1,...,N and cf> — $>/($>qN), being T the hopping matrix element between nearest neighbors. The time- 
dependent phase <f)t attached to each link accounts for the presence of the external magnetic flux. The term H a 
represents the lead and the reservoir. For the moment we do not give any explicit hamiltonian for this system. 
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It is enough to assume that it is an infinite system of non-interacting electrons with a chemical potential p. The 
hamiltonian describing the contact is 

Hcontact = -T la (c\c a + C^Cl). (3) 

In terms of Green functions, the current through the contact is written as 

J la = ^Re[T la G< a (t,t)}, (4) 

where Gf a (t,t') = i(c\(t)c a (t')). The latter obeys the following equation, 

G< (t,t') = -T la J dttiG^t^iti-t^ + Gf^t^ih-t')}. (5) 

The functions {t\ — t'),g^(ti — t') and g< (ti — i'), with g^(t — t') = [g^(t' — t)]* are equilibrium Green functions 
corresponding to the noninteracting system a before the coupling to the ring is performed. In terms of the spectral 
density p a (oj) = — 2Im[g^(oj)] and the Fermi function = 0(p — uj) of this system, they read 

g K(t-t') = -iG(t-t>) J |£p a ( w ) e -Mt-0 

~ t ')= i J ^U^P^e-^ 1 -^. (6) 

The Green functions G^ x (t, t') and GJ ; 1 (t,t') correspond to the ring coupled to the system a. They satisfy the 
following equations 



Gfi(M') = Si!i(M')+ / dfcidt 3 G^ 1 (t,t 1 )Ef(ti,-t 2 )^ 1 (t 2 ,t / ), 

G< A (t,if) = J dhdtzG^fotjEfih - t 2 )Gf A (t 2 ,if), (7) 

where g^i(t,t') is the retarded Green function for the uncoupled ring penetrated by the magnetic field. The self- 
energies 

^{t-t>) = \T la \ 2 g*(t-t') 1 

Z<(t-t>) = \T la \ 2 g<(t-t>) (8) 



account for the corrections due to the escape to the leads |1J,|13|] . This set of equations provide the exact solution for 
the time-dependent current J\ a (t). 

We now introduce one simple model to describe the system a. To simplify the numerical procedure, we assume that 
the ensuing electronic band has a large bandwidth and a constant density of states p a (u) = T. Thus, g£(t — t') ~ 
— iTS(t — t') and the following expression for the time-dependent current is obtained 

Jia(t) = f a J ^f a (u)(Im[G^(t,u)} + a\Gf tl (t,u)\ 2 ), (9) 

where a = |T lct | 2 r, with £^(t — t') = iaS(t — t'), while the Fourier transform for the retarded Green function is defined 



G« (*,«)= / diV^H'-^G* (*,*'), (10) 



as, 



with 77 = + . 

We remark now some salient issues related with expression ([)[). The function pi(t,u>) = — 2Im[Gf" l {t, ui)\ defines 
a generalized time-dependent spectral function for site I, while the function Ni(t,ui) = 2a\Gf 1 (t, ui)\ 2 defines the 
occupation of this site at time t, since 

Nt(t) = -iG<(t,t) = I ^/oWM), (11) 
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where Ni(t) = (clcj) is the mean number of particles at site I. In an equilibrium system or in a steady regime 
where G^it, u) does not have an explicit time dependence, it is easy to prove that Ni(ui) — pi{u). However, in a 
non-equilibrium time-dependent regime the local Green function Gfi(i,u;) depends explicitly on time and can develop 
a complex structure as a function of u>. Therefore, the two functions are in general different. This property is crucial 
for the physical effects we aim to address. In particular, the fact that in a stationary situation these two functions 
coincide, cause the current J\ a {t) to vanish. Instead, in a non-equilibrium regime, an ac current can flow through 
the lead. As we shall see, this current exhibits a periodic behavior which is originated by the existence of Bloch 
oscillations in the ring. 

To complete the theoretical description of the system, we define the current flowing through the bond between sites 
I and I + 1 along the ring 

Ju +1 = ^Re[Te-^G< +1 (t,t)], (12) 

where 

G< +1 (t,i) =i J ^aG^{t,uj)[G^ ltl {t^)]\ (13) 

with 

G*, n (t,t')=g* tn (t,t')+i J dhGK^tJagKjhJ). (14) 

From the practical point of view, the problem reduces to the calculation of the retarded Green functions n (t, t') 
and G^ n (t 7 t') corresponding the loop penetrated by the magnetic field unconected and coupled to the external 
reservoir, respectively. This is the subject of Sec. III. 



III. CALCULATION OF THE RETARDED GREEN FUNCTIONS 
A. The metallic ring threaded by a flux $m (£) = <&t 

We are interested in calculating the retarded Green function g^ n (t,t'), where m,n label the m-th and n-th sites 
of the ring, respectively. The usual procedure to compute this quantity is to take advantage of the translational 
invariance of the system 0. In this subsection, we shall discuss an alternative way to obtain exactly g^ n (t,t'). As 
we shall see, this route leads to a more transparent physical picture for the relevant energy scales of the problem while 
it has the additional advantage of being adequate to solve more general problems without translational invariance. We 
shall see in the next subsection that this procedure also provides a very convenient framework to treat the coupling 
of the ring with the external lead and reservoir. 

The Dyson equation for the retarded Green function leads to the following set of coupled differential equations 

~ f ) + Te**9%,n+i(t> + re- 4 **^* n -i(*. = *(* - W- ( 15 ) 

Our strategy is to perform a gauge transformation in order to translate the original problem into an equivalent one 
where the explicit time dependence appears only in the matrix elements involving the link between the first and the 
last sites of the ring. By defining the gauge transformation 

c„ = exp[irt0i]c n , (16) 

the Green function transforms as 

= exp[icj>(mt - nt')]g^ n (t, t'). (17) 
From (|l5|), it is found that g^ n (t,t') satisfies 
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-iKQp9*,n(t,t') - Hng^ n (t,t')+Tg^ n+1 (t,t') +Tg^ n _ 1 (t,t') 
= 5{t,t')S m , n , n?l,N 

-ih^g* A (t,t') - h4>g* A (t,t') +Tg^ 2 (t,t') +T N1 (t')g^ N (t,t') 
= S(t,t')S mA , n = l 

-^9^,nM - h §^9*At,t') +T 1N (t')gK A (t,t') +Tg* tN _ x {t,t') 

= S(t,t')S mtN . n = N, (18) 
where the only matrix elements that exhibit an explicit time dependence are 

T 1N (t) = exp[i— t]T, (19) 

with [T N1 (t)}* =T 1N (t). 

The above set can be interpreted as the Dyson equation for the retarded Green function corresponding to the 
problem of a chain with N sites of tight-binding electrons in the presence of a dc voltage 

Vi = M(f>. (20) 

The chain is closed at its ends by a time dependent link Tuy(t). We formally treat this system as indicated in the 
scheme of Fig. 2, by separating the ensuing hamiltonian in an "unperturbed" part that describes an open tight-binding 
chain, 

N-l N 

H° = fcfc+i + h.c) + Vicjci, (21) 

i=i i=i 

and a time-dependent "perturbation" 

Hwit) = -T 1N (t)c\c N - T N1 (t)4,ci, (22) 

which describes a time-dependent hopping between the first and the last sites of the chain. Within this picture, the 
integral form for the Dyson equation reads 



=0m,n(*.*O- / dt^ N {t, tl)T N1 {t X )gl n (t U t') 
Jt> 

dt 1 g^ 1 (t,t 1 )T 1N (t 1 )g° N Jt u t r ), (23) 



where the matrix elements Tjsn(t) and T\iq{t) define time-dependent instantaneous self-energies. The "unperturbed" 
Green function g^^t^t') = g„„(i — t') is an equilibrium Green function corresponding to the hamiltonian Hq. 
Although there is no trivial analytical expression for g^ n (t — t'), it can be easily obtained in a finite-size system from 
the numerical solution of hamiltonian (|2l]). In terms of the eigenvalues E v and eigenvectors \v) — J^i ^-flO °f -^Oj h 
can be written as 

JV 

<«(* - t') = -*0(t - f ) KiK exp[-^(i - t% (24) 

After performing in Eqs. ( ^3|) the Fourier transform defined in (|Io|), we obtain for n = 1 

i (t , w) + g*, N (t ,u + §- )T m (t)gl t (u) + g^ (t, lj - — )T X N (t)g% ;1 (w) = ^ (w) 
g*, N (t,cj) +g% N (t,w+ §-)T N1 (t)gl N (u) +^,i(t,w- s-JTijvW^H = ff ^(w). (25) 

For each time i, the solution of the above set of linear equations provides the complete exact solution of the problem. 
The Green function for arbitrary to, n is obtained from 
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g£, n (t,w) = <„(t,w) -3m,jv(*»^+ i-)Tm(t)gl n (u) -g* a (t,u>- —)T 1N (t)g° N!n (uj), (26) 
while the Fourier transform of the retarded Green function of the original problem is 

ffm,n(*, w ) = exp[i<^(m - n)t]g^ (t, uj - mf). (27) 



B. The ring coupled to the particle reservoir 

We now develop a computational strategy to obtain the retarded Green function G^ n (t, t 1 ) for the coupled system. 
The differential form of the Dyson equation (Oh for this Green function reads 



-^^G« n (t,0-w*n 1 iG» n + Te-** t '(?« wfl (t J * / ) 

+Te^ t 'G^ n _ 1 (t, t') = S m , n S(t - t'). (28) 

As in the previous subsection, we perform the gauge transformation ( |l6| ) in order to obtain a set of equations with 
a reduced number of time-dependent matrix elements. The Green function transforms as 

G* „(i, = exp[^(mt - nt')]G* f') (29) 
and G^ n (t,t') satisfies the following equations 

~^Gl. n {t,t') - hctmGl tn {t,t') +TGl tn+l {t,t') +T^ n _ 1 (t,t') 
= 5 m , n S(t-t') n^l,N 

-ih-^G* A (t,t') - HGl^t') - iaG* tl (t,t') +TG* t2 (t,t') +T N , 1 (t')G^ N (t,t') 
= 5 m ,i6(t - 1') 

-ih-^G^ N (t,t') - h^-G* tN (t,t') + T 1 Mt')G*, 1 (t,t') +TG^ jN _ 1 (t,t') 

= S m . N S(t - t') (30) 

The above set is formally identical to the Dyson equation of the problem illustrated in Fig. 3. This consists in an 
open tight-binding chain of N — 1 sites labeled by I = 2, . . . , N with potentials Vj (see (|20|)), plus an additional site 
labeled by / = 1 (indicated with a square in Fig. 3). The latter represents a site of the ring coupled to the lead and 
reservoir. It has a local potential V\ and a self-energy correction ia due to the coupling to the system a. This site is 
connected to one end of the chain through a hopping T and to the other end, through the time-dependent hopping 
Tijy(t) (|l9|). The integral form of the Dyson equation results, 

G* n (M') = g° m , n (t,t') - J dudt^Gt^tM^Ni^gl^t') - 

dwdhG^it^TiNit^glAh^-T J dtiE^t, tO - 

T J dtiG^t, *!)<??, n (ti,*0, (31) 
which by recourse to the Fourier transform ([To]) leads to the following set 

G^ 1 (t,cj)=g^ tl (u J )-Te- i ^ t G^ N (t,u J + n§-)gl^ 
G^ N (t,cj) =g° m<N (u J )-Te i ^ t G* A (t,u>-h§-)g% tN (u>)-^ 

Gl, 2 {t^) =< 2 M-Te i ^ t G^ 1 (t, W -^|-)^ 2 M-TG^ 1 (t )W ) 52 )2 M, (32) 
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where n (u>) with m, n = 2, . . . , N is the Green function of the chain formed with sites 2, . . . , N with a bias Vi, 
while g® 1 = 1/(uj — <p + iu) is the Green function for the remaining site in interaction with system a. 

After some rather simple algebra, it is easy to see that Eqs. ( p?^ ) can be reduced to the set (^5|) by replacing in the 
latter g*, n (t,u) by G m n {t,uj) and g° ln (t,uj) by G° m n (t,uj), being 

G&.nM = <„M - TG^o,)^ - TGl i2 ( W )4( W ). (33) 

A scheme of the underlying physical processes related with this equation is given in Fig. 3b. It corresponds to the 
coupling of the open chain with the site I — 1 by means of a single hopping element T. Equations ( P5|) correspond to 
the closing of the resulting A-site chain by means of a time-dependent hopping process (^) between its ends. 

In conclusion, the exact Green function for the loop penetrated by the time-dependent magnetic field, as an 
isolated system or as a system coupled to an external particle reservoir, is given by the solution of Eqs. (p5|), provided 
that a suitable expression for the "unpperturbed" Green function, g^ n (t,u>) or G^ n (f, u/) correspondingly, is given. 
We must emphasize that although in both cases the final set of equations is formally identical, physical processes 
of different nature are enclosed in the "unperturbed" Green function. In particular in the case considered in this 
subsection, the effect of the lead represented by the self-energy ia introduces inelastic scattering events that produce 
loss of coherence in the propagation of the wave function through the ring and dissipation. 

The imaginary part of the local Green function defines a generalized timc-dcpcndcnt density of states pi(t,u) — 
— 2Im[Gf' l {t, u>)], in terms of which the local Green function can be written as 

G^)=[^^^- (34) 
J Ztt oj — OJ -\- IT] 

This quantity is in general expected to be a periodic function on time, as suggested by the structure of the set (pH). 
However, we shall later see that in some situations it can also become independent of time. 



IV. RESULTS 



A. Dynamical regimes 

We first consider the case introduced in Section IIIA, which corresponds to the ring under the influence of the 
time-varying magnetic flux, but disconnected from the external wire and particle reservoir. Our aim is to understand 
first the dynamical behavior of this system, free from the effects of dissipation introduced by the coupling to the 
external reservoir. 

An important issue is that if g® N (u>) = g 1 ^ 1 (w) = 0, the equations in J25| ) can be decoupled. In that case, the 
solution becomes stationary, i.e. does not have any dependence on time and reads 

The latter property implies that the local Green functions of the original problem gi,i{t,u>) become independent of 
time, while the non-local ones exhibit a trivial time dependence introduced by the gauge transformation (|27|). 

In order to find out the physical conditions in which this stationary solution is expected, we study the behavior 
of g® n (lu). The relevant parameters are the magnitude HQ/Qq relative to the bandwidth of the tight-binding chain 
W = AT and the size of the system N. We numerically diagonalized the hamiltonian (|2l]) for different ratios h§/$>oW 
and studied systems with sizes up to N = 800 sites. Typical curves for the imaginary part 

N 

Im[gl N {u)] =-ttJ2 AtA» N S(u - E v ). (36) 

for a chain with N = 20 sites are shown in Fig. 4. It is clear that for all u, the magnitude of this quantity is 
significantly larger in the cases where TlQ/Qq < W. 

To be able to perform a more quantitative analysis of the behavior of g® N (w) as a function of N, we consider 
the function g m ax{N), which is defined as the maximum value achieved by |/m[g? i (a;)]| for a given size N. Plots of 
gmax(N) as a function of N reveal that there is a change in the scaling of this quantity as fr$/$oW varies. The case 
7t$ < <&$W is illustrated in Fig. 5. We found that the scaling law is 
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g max (N) tx ft$/$ < W. (37) 

The data of Fig. 5 can be adequately fitted with 7 ~ 1 for h$>/$ W = 0, 0.25, 0.5, 0.75 and 7 - 1.4 for h<f>/<S> W = 1. 
The opposite situation, where 7z.<E» / <E»o > is shown in Fig. 6. We found that in this case, the scaling law is 

N 

9ma X (N) oc e~— , 7i$/$ >W. (38) 

By fitting the numerical data it is obtained L a ~ 8 for ft^/'&oW^ = 1-25 and L Q ~ 4 for 7i<I>/$oW = 1.5. 

The existence of two kinds of behaviors in the scaling of g® n (uj) has a rather simple explanation. We recall that 
this Green function corresponds to a system of electrons in an open tight-binding chain with a bias 7j$ / $0 between 
its edges. The different scaling laws result from a concurrence between the kinetic energy of the electrons in the chain 
and the strength of the bias. In fact, the Green function g® N (t — t') is by definition related with overlaps of the 
form (if!i(t)\i$!' L {t')), where (*)) and |^((t)) are the wave functions of the system when an electron is created and 
destroyed, respectively, in site I at time t. This quantity is thus a measure for the ability of an electron to move from 
the first to the last site of the chain in a time t — t'. In the free chain, the probability of such a movement decreases as 
1/N as the size of the system increases, independently of how large the time it takes is assumed to be. Our numerical 
analysis indicates that this is also the case in the presence of a bias, provided that the strength of the bias is not 
larger than the bandwidth associated with the kinetic energy of the electrons. As the bias increases and overcomes 
the value of the bandwidth, this probability becomes smaller and decreases exponentially with the size of the system. 

Therefore, for magnetic fields varying at a rate such that h$/$>o < W, the Green functions 5^ n (^' w ) tend to the 
steady solution (|35| ) as TV -7 when the size of the system increases. Hence, for finite size systems, time dependent 
features are clearly observable. As a function of time, the functions g^ n (t,oj) oscillate with the period r = 2n& /$>. 
The latter being the period of the Bloch oscillations due to the induced electric field of magnitude F = $>/cL. Instead, 
for rapid variations of the magnetic field such that h$/$o > W, the Green functions tend to ( |3^ ) exponentially fast 
as the length of the circumference increases. As a consequence, within this regime, time dependent features can be 
observed only in very small rings of lengths smaller than L a . 

We close this subsection with some comments about the behavior introduced by inelastic scattering processes in the 
dynamical behavior above mentioned. We have shown in section IIIB that the Green function for the coupled ring is 
also given by the solution of Eqs. ( p5| ) with the replacement of (u) by the function G^ n (w) which contains the 
corrections due to the escape to the leads. Thus, the equations in d2q) are decoupled and a stationary solution 

— R — R ^mn(^) 

G m ,n (W) = G m ,n M = l_ T 2 G O i(w)G O rN(a; + $/$o) ■ ( 39 ) 

exists also in this case when the function G\ L {u) becomes vanishingly small. As before, this quantity is related with 
the probability of an electron to travel from the first to the last sites of the chain. In this case, the new ingredient is 
the existence of a typical inelastic scattering length Li n beyond which the coherence in the wave packet propagation 
is loss. The latter is related to the imaginary part of the self-energy through a oc l/Li n which generates damping 



and produces an exponential decay ~ exp[— aL/Li n ] in the Green functions G^ n (w) |15|. Therefore, for rings with 
circumferences L < Li n the same two dynamical regimes defined for the disconnected ring can be distinguished. 
Instead, as soon as L > Lj„, inelastic scattering due to the presence of the external reservoir dominate the physical 
behavior of the system, the Green functions coincide with the steady solution ([35]) and no time dependent features 
are observed. An example of significant time-dependent behavior is shown in Fig. 7 where typical shapes of the local 
generalized density of states px(t, ui) for a ring of N = 20 sites with a = 0.025W at different times are depicted. As a 
function of u>, p\(t,uj) can be either positive or negative as a consequence of the far from equilibrium nature of this 
regime. The sum rule 

— f>i(t,u) = l (40) 

Z7T 

is, however, satisfied, justifying the name of density of states. 



B. The current through the lead 

Although a classical textbook analysis indicates that the current through the single lead coupling the ring to the 
particle reservoir should vanish, within the time-dependent regime defined in the previous subsection, there is an ac 
contribution which could be measured in a real experiment. 
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The time-dependent regime is characterized by an imaginary part of the self energy a leading to an inelastic 
scattering length Li n larger than the size of the ring and a variation of the magnetic field satisfying the condition 
UQ/QoleqW. The propagation of the wave packet maintains the coherence, the Green functions gi,i(t, u) and C?i,i(t, u) 
oscillate with the period r = 27r$o /$ of Bloch oscillations driven by the induced electric field, and so does the current 
through the lead. Examples for a ring with N = 20 sites are shown in Fig. 8 where Ji a (t), obtained from (|^), is shown 
for two different values of <I>. The imaginary part of the self-energy is taken to be a — 0.2VF. The physical picture 
that emerges is that the induced electric field produces dramatic changes in the generalized time-dependent density 
of states in the ring p\[t, u>). This quantity is periodic in time and eventhough it satisfies the sum rule (|40|), it can 
achieve positive as well as negative values as a function of u> (see Fig. 7). In particular, this function does not satisfy 
the property of being equal to the occupation function N\ (i, u>) . As a consequence the integral (^) does not vanish in 
general but displays oscillations with zero mean value as a function of time. A somehow similar phenomenon in the 
context of optical properties of semiconductors is the so called Franz-Keldysh effect JL6| |. In that case, the changes in 
the density of states due to the influence of time-dependent fields cause the absorption coefficient to become finite for 
photon energies bellow the band-edge. 

Interestingly, due to the symmetry of the spectrum (see Fig. 7) the current J\ a {t) is an antisymmetric function 
of the chemical potential \x. This property is due to the particular geometry of the system we are considering. Any 
spectral asymmetry, introduced, for example, by the presence of barriers and/or disorder in the ring should remove 
this feature. 

The effect of inelastic scattering is illustrated in Fig. 9, where J\ a is shown for several values of a. Inelastic 
scattering introduces damping and Li n decreases as a increases. In terms of the spectrum, the peaks of pi(t, u>) get 
wider and, for very large a, the time dependent features tend to be wiped out. On the other hand, very weak inelastic 
scattering, which occurs for example for a very small value of Ti Q , tends to forbid the quantum tunneling between 
the ring and the reservoir. Hence, the amplitude for the oscillations of J\ a is vanishingly small for small a, tends to 
decrease for very large a and attains a maximum in between. 



C. The current along the ring 

The resistive effect caused by inelastic scattering in metallic loops was studied some time ago by Landauer and 
Biittiker ||, who proposed a phenomenological equation of motion to analyze the relaxation processes in a ring 
penetrated by a time-dependent magnetic flux. A dc component in the electric current along the ring was found to 
increase as inelastic scattering is introduced and to decrease as inelastic scattering becomes sufficiently intense. This 
analysis was later supplemented by Biittiker |[(|, who proposed a single lead connecting the ring with an external 
reservoir like in Fig. 1, as a simple device to generate a resistance in the metallic loop. 

In this section we show that the theoretical approach adopted in the present work is also appropriate to investigate 
the effect of inelastic scattering in the current along the ring. Results for the current (|l^) through the bond between 
the sites I = 1,2 of the ring are shown in Fig. 10 for a system with N = 20 sites. We consider parameters within 
the time-dependent regime. The imaginary part of the self-energy is assumed to be a — 0.2VF. As a function of 
time, the current Jn(t) displays oscillations with the period r = 27r<E>o / of Bloch oscillations. Note that the factor 
e -i</>t m (12) exactly cancels the time-dependent exponential introduced by the gauge transformation (p7]). Thus, 



the oscillations of the current Jz,;+i result as a consequence of the time-dependent behavior of the Green functions 

(*m n (t, lo) obtained from the solution of Eqs. ( |25| ) and disappear in the steady regime. The remarkable feature is the 
existence of a dc component in the current, which depends on the chemical potential of the reservoir, as well as on 
the rate ?i$/$o for fixed a. The effect of a on the magnitude of the dc component is shown in Fig. 11, where it is 
observed the behavior predicted in Ref namely the magnitude of this component is vanishes for small as well as 
for very large a and has a maximum for intermediate strength of dissipation. 

To close this section, we show the behavior of the mean occupation number Ni(t) along the ring. The latter is 
calculated from ([ll]) for a ring with TV = 20 sites. Results are shown in Fig. 12 for some values of I. As a function of 
time it displays oscillations with period r around the mean value 

(Ni) = i J dtNi{t), (41) 

that depends strongly on the position /, as shown in Fig. 13. The behavior of (JVj) is a measure of the effective mean 
potential along the chain. The latter exhibits a significant drop near the site / = 1 where the lead is connected, which 
becomes more abrupt as the strength of the inelastic scattering introduced by a increases. 
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V. SUMMARY AND CONCLUSIONS 



We have studied the transport properties of a ring threaded by a magnetic flux with a linear dependence on time 
^Ai(t) = <£>£. We considered an experimental setup to investigate the transport properties of this system, which consists 
in a single lead connecting the ring to an external particle reservoir. We employed a theoretical approach based in 
Baym-Kadanoff-Kcldysh non-equilibrium Green functions. We derived a system of linear equations to calculate in 
real space the retarded Green function of the disconnected ring as well as for the ring coupled to the particle reservoir. 

This method provides a transparent picture for the different energy scales governing the dynamics of the system. 
In particular, we identified two kind of behaviors in the Green functions of the uncoupled ring, depending on the ratio 
K<f>/<f>oW. For magnetic fields varying at a rate such that ?i<I>/$o < W , the time dependent features are important 
in the behavior of the Green functions of finite rings and tend to disappear as iV~ 7 when the size of the system 
increases. Instead, for ?i$/$o > W, the Green functions tend exponentially fast to a static solution as the size of 
the system increases. These two kinds of behaviors manifest themselves in the transport properties of the system 
when the inelastic scattering length is larger than the perimeter of the ring. Within the time-dependent regime, the 
current through the lead as well as the current along the ring and the charge inside it, display oscillations with the 
period r = 27r<&o / ^ of Bloch oscillations driven by the induced electric field. In the situation of a rapid magnetic field 
variation satisfying ft$/$o > W, no current is measured in the lead and no oscillatory behavior is observed except 
for extremely small rings. We stress that this effect has its grounds in the concurrence between the kinetic energy 
of the electrons and the strength of the dc voltage due to the induced electric field, being independent of the loss of 
coherence introduced by inelastic scattering. 

We finally analyzed the effect of inelastic scattering caused by the coupling with the external reservoir in the time- 
dependent transport. We found that the amplitude of the oscillations of the current through the lead as well as the 
dc component of the current along the ring, tend to vanish in the two limiting situations of very weak and very strong 
inelastic scattering, attaining a maximum for intermediate strengths. 

On the technical side, we find that the present approach provides an useful theoretical tool to investigate the 
effects of time dependent magnetic fluxes and induced electric fields in other interesting systems of larger complexity, 
like metallic loops with disorder, rings with inserted barriers and quantum dots and one dimensional systems with 
electronic correlations. 
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FIG. 1. (a) Scheme of the experimental setup to study the charge transport in a metallic ring threaded by a time-dependent 
magnetic field, (b) Theoretical model for the device in (a) . 



T 1N (t) 




FIG. 2. Scheme of the physical processes involved in the Dyson equation for the retarded Green function (C 



[17] The analytical expresion for the Fourier transform is g^(t,t') — —i®(t — t') J* ds exp[— iek(s)], where £fc(s) = — 2Tcos(fc + 
<f>s), with k = 2mv/N, n = 0, . . . , N - 1. 
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(a) 




FIG. 3. Scheme of the physical processes involved in the Dyson equation for the retarded Green function. Figs, (a) and (b) 
correspond to Eqs. (B3) and (B3I), respectively. 
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FIG. 4. Imaginary part of the Green function gi N (co) for a system of N = 20 sites and magnetic fluxes varying at rates 
7i$/$o = 0,0.75W, 1.25W, 1.5W for (a), (b), (c) and (d), respectively. 
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as a function of 1/N. Circles, squares, diamonds and triangles correspond to h<fr/&o = 0.25W, 0.33W, 0.5W, W, 





FIG. 7. The generalized time-dependent density of states for a chain of N = 20 sites with = 0.5W a — 0.025W 7 at 

different times t = 0, 1, 2 (a, b and c, respectively). 
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FIG. 8. The current Ji a (t) between the ring and the reservoir as a function of time for a ring with L = 20 sites and a 
self-energy a — 0.2W. Circles correspond to S,4>/«3>o = 0.25W for a chemical potential of the reservoir fi — 0.25W (solid circles) 
and fj, — — 0.25W (open circles). The solid line corresponds to ?z<&/ <E>o = 0.5W with fi = — 0.25W. 
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FIG. 9. The current Ji a (t) between the ring and the reservoir as a function of time for a ring with L — 20 sites. The 
chemical potential of the reservoir is fi — —0.25W and ?j,4> / <3>o = The thin solid, dashed, dot-dashed and thick solid lines, 

filled circles, opaque circles and filled squares correspond to a — 0.0025W, 0.025W, 0.2W, 0.5H 7 , W, 2.5W, 12.5VK, respectively. 
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FIG. 11. The current Ji2(t) along a ring with L — 20 sites and ?i$/$o = 0.5W and /i = 
squares, open circles and open squares correspond to a = 0.0025W, 0.2W, W, 2.5W, respectively. 
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